Possible Anderson transition below two dimensions in disordered systems of 

noninteracting electrons 
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We investigate the possibility of an Anderson transition below two dimensions in disordered 
systems of non-interacting electrons with symplectic symmetry. Numerical analysis of energy level 
statistics and conductance statistics on Sierpinski carpets with spin-orbit coupling indicates the 
occurrence of an Anderson transition below two dimensions. 



According to the scaling theory of Abrahams et al. for 
non-interacting electrons in disordered systems, all states 
are localized in two dimensions QD) and the Anderson 
transition occurs only above 2D This prediction ap- 
plies to systems with orthogonal symmetry, i.e., systems 
with time reversal and spin rotation symmetry [H |^ Q . 
Accordingly, it is believed that the lower critical dimen- 
sion for the Anderson transition in systems with orthog- 
onal symmetry is 



,(orth) 



(1) 



The prediction of Abrahams et al. does not apply to 
systems with symplectic symmetry, i.e. in systems with 
time reversal symmetry but in which spin rotation sym- 
metry is broken by spin-orbit coupling 0, 0| ■ For such 
systems it is known that there is a transition in 2D, if 
the spin-orbit interaction is sufficiently strong. 

The prediction of Abrahams et al. rests on an argu- 
ment concerning the form of the /3 function that describes 
the scaling of the conductance. Reasonable assumptions 
concerning the asymptotic behavior of this function in 
the strongly metallic and localized limits, respectively, 
and the assumption that the /3 function is monotonic lead 
to their conclusion. As we explain below it is this latter 
assumption that does not hold in systems with symplec- 
tic symmetry. 

In Fig. n] we show the B function that describes the 
scaling of the quantity AjSl 



/3(lnA) 



dlnA 
dlnL' 
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Here A is the ratio of the quasi- ID localization length to 
the system width for electrons on a long quasi- ID system 
of width L. The (3 functions that describe the scaling of 
different physical quantities are expected to differ in de- 
tail, but for systems with the same dimensionality and 
symmetry all are expected to share some common prop- 
erties. For example, a zero of the (3 function indicates a 




FIG. 1: The numerically estimated (3 functions for systems 
with symplectic symmetry in 2D and 3D (solid lines) |^ 
as well as the schematic ones for systems with orthogonal 
symmetry (dotted lines). 



transition and this property should be common to all the 
/3 functions. Also the slope at the zero, which is related 
to the critical exponent v describing the divergence of 
the localization length, should be the same for all the j3 
functions. 

The solid lines in Fig. ^ are numerical estimates of 
the function (jSJ for systems with symplectic symmetry 
0, Q- The (3 function tends to d — 2, with d being the 
dimensionality, in the metallic limit and becomes nega- 
tive in the localized limit. In contrast to systems with 
orthogonal symmetry (shown schematically, see Ref. 
for numerical result) the (3 functions are non-monotonic 
and a fixed point exists even in 2D. It is thought that the 
13 function for the conductance behaves similarly Q . 

Looking at Fig. ^the obvious question is, "What hap- 
pens below 2D for systems with symplectic symmetry?" 
Does the (3 function change discontinuously and suddenly 
become monotonic? Or does it change smoothly? In the 



2 



latter case, we expect the Anderson transition to per- 
sist below 2D. This is a long standing problem Our 
numerical results suggest that the Anderson transition 
persists below 2D for systems with symplectic symme- 
try. It is known that states in (quasi-) ID systems with 
symplectic symmetry are always localized, so the lower 
critical dimension must be at least one. This would mean 
that 



1 < d 



(symp) 



< 2. 
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To study this problem we have simulated non- 
interacting electrons on the Sierpinski carpet with 
spin-orbit coupling. The Sierpinski carpet is constructed 
by repeating the following iteration. We start from abxb 
square lattice with a central c x c square removed. This 
first generation or "generator" is denoted as SC(&, c, 1). 
At each step the next generation is constructed by mag- 
nifying the lattice b times and replacing each site of the 
current generation with the generator. The linear size of 
the fcth generation SC(6, c, k) is L = b*' and the number 
of sites is = [b'^ — (?)^ . The Sierpinski carpet be- 
comes a true fractal in the mathematical sense when the 
generation number k tends to infinity. This true frac- 
tal is denoted as SC(6, c). According to a lot of stud- 
ies of thermal magnetic phase transitions, the Sierpin- 
ski carpets are useful to study phase transitions in frac- 
tal dimensionality less than two |0| (see also a recent 
work Ref. 0| and references therein). We have studied 
SC(3, 1, k) and SC(5, 1, k) with fixed boundary conditions 
imposed in both directions. 

A fractal system has a fractal dimension df. This di- 
mension describes how the mass of the fractal depends 
on the size of the fractal. The fractal dimension df of the 
Sierpinski carpet SC(6, c), defined by A^ = i"*' , is equal 
to df = ln(62 - c2)/ln6. We find d{ « 1.893 for SC(3, 1) 
and df w 1.975 for SC(5jJ_). A fractal system also has 
a spectral dimension ds |l4l |. This dimension describes 
the spectra of the low energy vibrations on the fractal. 
Numerical studies of the Anderson transition in bifractal 
systems indicate that this is the relevant dimension when 
discussing the Anderson transition 0, . The spectral 
dimension of the Sierpinski carpet has been estimated 
by simulating a classical random walk on the fractal . 
The values reported for SC(3, 1) and SC(5, 1) are, respec- 
tively, 4 = 1-785 ± 0.008 and 4 = 1.940 ± 0.009. 

We have simulated the SU(2) model [illll. 
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Here c\^{cia) denotes the creation (annihilation) oper- 
ator of an electron at the site i with spin a. This 
model has symplectic symmetry. There is an on-site 
random potential which is identically and indepen- 
dently distributed with uniform probability in the range 
[—W/2,W/2]. There is also nearest neighbor hopping 



between lattice sites with the hopping matrices R{i,j) 
sampled uniformly from the group SU(2). Apart from 
the fractal lattice the model is exactly as described in 
0| and the reader should refer to these for further 
details. 

To investigate the localization of electrons on these 
fractals we have examined the energy level statistics of 
Q). In particular, we have looked at the distribution P{s) 
of the nearest neighbor level spacing s measured in units 
of the mean level spacing It is known that P(s) is 
well approximated by the Wigner surmise for the Gaus- 
sian symplectic ensemble -Pgse(s) = As'^e^-^'' in the ex- 
tended limit, and the Poisson distribution Pp{s) = e~^ 
in the localized limit. We have calculated the energy 
eigenvalues using either LAPACK or the Lanczos method 
[igj . An unfolding procedure has been applied to the level 
spacings to compensate for the variation with energy of 
the ensemble average density of states. For a quantitative 
analysis of P(s), we define 



J^" dsPosEis) - J^" dsPpisY 
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In the extended limit Y^^ — 1, and in the localized limit 
Ys„ — 0. Note that since each energy level has a twofold 
Kramers' degeneracy we consider only distinct eigenval- 
ues when calculating the spacing distribution. 

We have simulated lattices SC(5,l,fc) up to fc = 
4. Calculation of higher generations is not possible at 
present because the CPU time required is too large. For 
the first generation k = 1 the number of eigenvalues is 
small. Therefore, we restricted our analysis to a single 
pair of consecutive levels per sample, the smallest positive 
eigenvalue and the largest negative eigenvalue. When we 
choose the pair of levels in this special way, the probabil- 
ity P{s) of a level spacing with a value s is equal to sP{s). 
Therefore we multiply P{s) by 1/s to obtain level spac- 
ing distribution P{s). We simulated 1000000 samples for 
k = 1. For higher generations the number of eigenvalues 
is greater so we reverted to the usual procedure of ana- 
lyzing a sequence of consecutive levels. The number of 
samples was reduced accordingly. We used levels in the 
interval E = [-0.3,0.3]. We simulated 15000, 600 and 
10 samples for k =2, 3 and 4, respectively. The total 
number of the level spacings for each of the (fc, W) pair 
is thus of order 10^ to 10^ 

Figure |2 shows Y1.5 as a function of disorder W for 
SC(5, 1, k). There are large corrections to scaling in the 
first generation. The data for the higher generations fc = 
2, 3 and 4 indicate that an Anderson transition occurs on 
SC(5, 1) around W ~ 5. When disorder is weak enough, 
Y1.5 increases with fc for k > 2 indicating a delocalized 
phase. When disorder is strong enough, Y1.5 decreases as 
fc increases indicating a localized phase. 

Further analysis of the level statistics on SC(5,1) is 
based on the assumption that finite size scaling applies 
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TABLE I; The details of the scahng analyses for SC(5, l,fc). Data with W in [4.2,5.7] have been used. The range of Ysg 
has also been restricted as indicated. Here, A'p is the number of parameters, A^d the number of data, Q the goodness of fit 
probability, and Yc the critical value of Yag. The precision of the fitted parameters is expressed as 95% confidence intervals. 
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FIG. 2: Fi.s vs. W for SC(5, l,fc) (fc = 1,2,3,4). The solid 
lines are the best fits of (|7J to the data inside the box (shown 
by dashed lines). Outside the box the fit deviates from the 
numerical data because the expansion parameters il)!}!^ and 
w are no longer small. 

to fractals in the same way as to systems with integral 
dimension. Accordingly the statistic Yg^ should, in the 
absence of any corrections to scaling, obey the single pa- 
rameter scaling law 

n„ - /± - ^0 (^i'/"-) . (6) 

Here ^ is the localization(correlation) length, v the crit- 
ical exponent for the divergence of ^, and i\) a smooth 
function of disorder W that crosses zero linearly at the 
critical disorder Wc- The subscript ± distinguishes the 
scaling function in the delocalized and localized phases. 
It is, however, clear from Fig. [21 that corrections to scal- 
ing are not negligible. Taking account of corrections up 
to the first order in an irrelevant variable (?5>L^ we have 

US 111113 

Ysg = Fo (Vii/"-) + ^Lyp, (^Li/-) . (7) 

For the purpose of fitting, the scaling functions Fq and 
Fi are approximated by a power series up to the order 
nn and n/ in ij}L^^'^ . The functions ip and are also 
approximated by expansions in terms of dimensionless 



FIG. 3; The numerical data Y1.5 after subtraction of the 
corrections to scaling plotted as function of the ratio i/^. 

disorder w — (Wc — W)/Wc up to the order mjf and 
mj. The best fit to data is determined by minimizing 
the statistics and the precision of the parameters are 
determined using the Monte Carlo method [l^ . 

The best fit is shown in Fig. El and the results of the 
scaling analyses for so = 0.5 and so = 1-5 are tabulated 
in Table m The estimates of the critical disorder Wc and 
the critical exponent v for so = 0.5 and sq = 1.5 are con- 
sistent as required. Our estimate 1/ = 3.4±0.2 for SC(5,1) 
is clearly different from the value ly — 2.746 ± 0.009 in 
2D 6], reflecting the difference of the dimensionality. To 
exhibit single parameter scaling the data are re-plotted, 
after subtraction of the corrections due to an irrelevant 
variable, in Fig. (31 The two branches correspond to the 
delocalized and localized phases. 

To complement our study of the level statistics, we 
have also analyzed statistics of the Landauer conduc- 
tance g (in units of e^/h) [21I E3 . i24j . Two perfect leads 
of width L = are attached to the Sierpinski carpet 
SC(5, 1, k) and the recursive Green's function method is 
used to calculate the conductance 25] . The hopping ma- 
trices in the direction transverse to the current are set 
to the unit matrix. As a scaling quantity, we choose the 
typical conductance 

<?typ=e<'"«>. (8) 
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FIG. 4: The typical conductance gtyp vs. W for SC(5,l,k) 
{k — 2, 3, 4). The hnes are a guide to the eye only. 

Here () means the ensemble average. We have set the 
Fermi energy to -E = and have accumulated 3000 sam- 
ples for k = 2,3 and from 150 to 300 samples for k — 4. 
Figure ^ shows the typical conductance gtyp as a func- 
tion of W for SC(5,l,fc). These data also indicate the 
occurrence of an Anderson transition at about W ~ 5. 

We have also analyzed energy level statistics and con- 
ductance statistics on SC(3, up to the fourth gen- 
eration. We did not find any evidence for a delocalized 
phase. We cannot conclude from this that SC(3, 1) is be- 
low the lower critical dimension because the SU(2) model 
has significant randomness even when the on-site random 
potential is zero. 

In Fig. I^lwe have plotted the available results for the 
critical exponent for systems with orthogonal and sym- 
plectic symmetries. For 2D and 3D system with sym- 
plectic symmetry the critical exponent are estimated to 
he u ^ 2.746 ± 0.009 in 2D [| and = 1.37 ± 0.02 in 
3D 0. For systems with orthogonal symmetry we have 
used the estimates reported in ^ . In the figure we 
plot l/iy as a function of ds- 

We recall that the slope of the f3 function at the fixed 
point is equal to l/v. Assuming that the /? function 
changes continuously as the dimensionality changes, then 
just above the lower critical dimension the fixed point 
should correspond to the maximum of the (3 function. 
We therefore expect the critical exponent i' to diverge 
at the lower critical dimension. This is consistent with 
the theoretical analysis of the nonlinear a model with 
orthogonal symmetry where it was found that = 1/e 
for e = d — 2 ^ 1 [2^. The numerical data for orthogonal 
systems shown in Fig. are consistent with this. 

For systems with symplectic symmetry fewer data are 
available 't'I . When only the estimates for v in 2D and 
3D are considered, it can leave the impression (cf. the 
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FIG. 5: 1/u vs. ds for systems with orthogonal and symplec- 
tic symmetries. The values for 2D and 3D symplectic systems 
are from 6] and 0, respectively. The estimates for systems 
with orthogonal symmetry with dimension 2 < da < 3 are 
from 4] and that in 3D is from "5^. (The exponent u for 
fractals with 2 < ds < 3 was also estimated in Ref. ^3). The 
solid line indicates u — l/{ds — 2). 

dotted line in Fig. 01 that the lower critical dimension is 
1 in this case. However, our estimate of the critical expo- 
nent on SC(5, 1) is well below the dotted line, suggesting 
that the lower critical dimension is closer to 2 than to 1. 

To confirm the occurrence of the Anderson transition 
below 2D for systems of non-interacting electrons in dis- 
ordered systems with symplectic symmetry, and to deter- 
mine more precisely the lower critical dimension, further 
numerical studies on a variety of fractals over a larger 
range of generations are needed. Analysis of fractals 
with ds intermediate between the values for SC(3,1) and 
SC(5,1) might be particularly helpful. The nature of the 
possible delocalized phase also warrants further study. 

Recently electron transport through systems with frac- 
tal perimeter was studied experimentally Electron 
transport through fractal structures can be an interesting 
topic for further theoretical and experimental investiga- 
tion. 
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